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Abstract The dynamics of partons, hadrons and strings in relativistic nucleus- 
nucleus collisions is analyzed within the novel Parton-Hadron-String Dynamics 
(PHSD) transport approach, which is based on a dynamical quasiparticle model 
for partons (DQPM) matched to reproduce recent lattice-QCD results - including 
the partonic equation of state - in thermodynamic equilibrium. The transition from 
partonic to hadronic degrees of freedom is described by covariant transition rates 
for the fusion of quark-antiquark pairs or three quarks (antiquarks), respectively, 
obeying flavor current-conservation, color neutrality as well as energy-momentum 
conservation. The PHSD approach is applied to nucleus-nucleus collisions from low 
SIS to RHIC energies. The traces of partonic interactions are found in particular in 
the elliptic flow of hadrons as well as in their transverse mass spectra. 
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1 Introduction 

The 'Big Bang' scenario implies that in the first micro-seconds of the universe the 
entire state has emerged from a partonic system of quarks, antiquarks and gluons 
- a quark-gluon plasma (QGP) - to color neutral hadronic matter consisting of in- 
teracting hadronic states (and resonances) in which the partonic degrees of freedom 
are confined. The nature of confinement and the dynamics of this phase transition 
has motivated a large community for several decades and is still an outstanding 
question of todays physics. Early concepts of the QGP were guided by the idea of 
a weakly interacting system of partons which might be described by perturbative 
QCD (pQCD). However, experimental observations at the Relativistic Heavy Ion 
Collider (RHIC) indicated that the new medium created in ultrarelativistic Au+Au 
collisions is interacting more strongly than hadronic matter and consequently this 
concept had to be severely questioned. Moreover, in line with theoretical studies 
in Refs. 0] [2] the medium showed phenomena of an almost perfect liquid of 
partons flUE] as extracted from the strong radial expansion and the scaling of ellip- 
tic flow V2{pr) of mesons and baryons with the number of constituent quarks and 
antiquarks J4). 

The question about the properties of this (nonperturbative) QGP liquid is dis- 
cussed controversially in the literature and dynamical concepts describing the for- 
mation of color neutral hadrons from colored partons are scarce. A fundamental 
issue for hadronization models is the conservation of 4-momentum as well as the en- 
tropy problem, because by fusion/coalescence of massless (or low constituent mass) 
partons to color neutral bound states of low invariant mass (e.g. pions) the number 
of degrees of freedom and thus the total entropy is reduced in the hadronization pro- 
cess. This problem - a violation of the second law of thermodynamics as well as the 
conservation of four-momentum and flavor currents - has been addressed in Ref. J6| 
on the basis of the DQPM employing covariant transition rates for the fusion of 
'massive' quarks and antiquarks to color neutral hadronic resonances or strings. In 
fact, the dynamical studies for an expanding partonic fireball in Ref. [Q suggest that 
the these problems have come to a practical solution. 

A consistent dynamical approach - valid also for strongly interacting systems - 
can be formulated on the basis of Kadanoff-Baym (KB) equations [7j] or off-shell 
transport equations in phase-space representation, respectively (7). In the KB the- 
ory the field quanta are described in terms of dressed propagators with complex 
selfenergies. Whereas the real part of the selfenergies can be related to mean-field 
potentials (of Lorentz scalar, vector or tensor type), the imaginary parts provide in- 
formation about the lifetime and/or reaction rates of time-like 'particles' f8]. Once 
the proper (complex) selfenergies of the degrees of freedom are known the time evo- 
lution of the system is fully governed by off-shell transport equations (as described 
in Refs. QUI)- The determination/extraction of complex selfenergies for the par- 
tonic degrees of freedom has been performed before in Ref. [0 by fitting lattice 
QCD (1QCD) 'data' within the Dynamical QuasiParticle Model (DQPM). In fact, 
the DQPM allows for a simple and transparent interpretation of lattice QCD results 
for thermodynamic quantities as well as correlators and leads to effective strongly 
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interacting partonic quasiparticles with broad spectral functions. For a review on 
off-shell transport theory and results from the DQPM in comparison to 1QCD we 
refer the reader to Ref. (SJ- 

The actual implementations in the PHSD transport approach have been presented 
in detail in Refs. ifTOlfTTl . Here we present results for transverse mass spectra and 
elliptic flow of hadrons for heavy-ion collisions at relativistic energies in comparison 
to data from the experimental collaborations. 



2 The PHSD approach 

The dynamics of partons, hadrons and strings in relativistic nucleus-nucleus colli- 
sions is analyzed here within the Parton-Hadron-String Dynamics approach |6| [TO] 
[TT1 . In this transport approach the partonic dynamics is based on Kadanoff-Baym 
equations for Green functions with self-energies from the Dynamical QuasiParti- 
cle Model (DQPM) [9:| which describes QCD properties in terms of 'resummed' 
single-particle Green functions. In Ref. |TT|, the actual three DQPM parameters for 
the temperature-dependent effective coupling were fitted to the recent lattice QCD 
results of Ref. 1121 . The latter lead to a critical temperature T c w 160 MeV which 
corresponds to a critical energy density of E c « 0.5 GeV/fm 3 . In PHSD the parton 
spectral functions pj (j = q,q,g) are no longer 8 functions in the invariant mass 
squared as in conventional cascade or transport models but depend on the parton 
mass and width parameters: 

separately for quarks/antiquarks and gluons (J = q,q,g). With the convention E 2 (p 2 ) = 
p 2 +M 2 — Yj, the parameters Af? and Jj are directly related to the real and imaginary 
parts of the retarded self-energy, e.g. Ilj = M 2 — 2ijj(0. The spectral function (Hi is 
antisymmetric in (0 and normalized as 

f°° da> . . [°° d(0 „ . , 
J -co 2n a p ''( ffl ' p ) = J 2n 2(a Pj( m >v) = 1 ■ (2) 

The actual parameters in Eq. (03, i.e. the gluon mass M g and width y g - employed 
as input in the PHSD calculations - as well as the quark mass M q and width y q , are 
depicted in Fig.[T]as a function of the scaled temperature T/T c . As mentioned above 
these values for the masses and widths have been fixed by fitting the lattice QCD 
results from Ref. lfT2l in thermodynamic equilibrium. 

One might worry whether the quasiparticle properties - fixed in thermal equilib- 
rium - also should be appropriate for out-off equilibrium configurations. This ques- 
tion is nontrivial and can only be answered by detailed model investigations e.g. 
on the basis of Kadanoff-Baym equations. We recall that such studies have been 
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Fig. 1 The effective gluon mass 
M g and witdh j„ as function 
of the scaled temperature T /T c 
(red lines). The blue lines show 
the corresponding quantities for 
quarks. 




summarized in Ref. JH) for strongly interacting scalar fields that initially are far 
off-equilibrium and simulate momentum distributions of colliding systems at high 
relative momentum. The results for the effective parameters M and y, which corre- 
spond to the time-dependent pole mass and width of the propagator, indicate that 
the quasiparticle properties - except for the very early off-equilibrium configuration 
- are close to the equilibrium mass and width even though the phase-space distribu- 
tion of the particles is far from equilibrium (cf. Figs. 8 to 10 in Ref. O). Accord- 
ingly, we will adopt the equilibrium quasiparticle properties also for phase-space 
configurations out of equilibrium as appearing in relativistic heavy-ion collisions. 
The reader has to keep in mind that this approximation is far from being arbitrary, 
however, not fully equivalent to the exact solution. 

We recall that the DQPM allows to extract a potential energy density V p from 
the space-like part of the energy-momentum tensor which can be tabulated e.g. as 
a function of the scalar parton density p v . Derivatives of V p with respect to p s then 
define a scalar mean-field potential U s (p s ) which enters the equation of motion for 
the dynamical partonic quasiparticles. As one can see from Fig. [2 the scalar poten- 
tial is rather large and nonlinearly increases with p s . This implies that the repulsive 
force due to U s (p s ) will change in a non-monotonous way with the scalar density. 
The vector mean-field potential is not negligible, too, especially at high p s and in- 
duces a Lorentz force for the partons. Note that the vector mean-field vanishes with 
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decreasing scalar density whereas the scalar mean-field approaches a constant value 
for p. s -K). 

Furthermore, a two-body interaction strength can be extracted from the DQPM as 
well from the quasiparticle width in line with Ref. [3]. The transition from partonic 
to hadronic d.o.f. (and vice versa) is described by covariant transition rates for the 
fusion of quark-antiquark pairs or three quarks (antiquarks), respectively, obeying 
flavor current-conservation, color neutrality as well as energy-momentum conserva- 
tion ifTOirfTl . Since the dynamical quarks and antiquarks become very massive close 
to the phase transition, the formed resonant prehadronic color-dipole states (qq or 
qqq) are of high invariant mass, too, and sequentially decay to the groundstate me- 
son and baryon octets increasing the total entropy. 

On the hadronic side PHSD includes explicitly the baryon octet and decouplet, 
the 0~ - and 1 ~ -meson nonets as well as selected higher resonances as in the Hadron- 
String-Dynamics (HSD) approach lfl3l [141 . The color-neutral objects of higher 
masses (>1.5 GeV in case of baryonic states and >1.3 GeV in case of mesonic 
states) are treated as 'strings' (color-dipoles) that decay to the known (low-mass) 
hadrons according to the JETSET algorithm [13] ■ We discard an explicit recapitula- 
tion of the string formation and decay and refer the reader to the original work |fT31l . 
Note that PHSD and HSD (without explicit partonic degrees-of-freedom) merge at 
low energy density, in particular below the critical energy density e c w 0.5 GeV/fm 3 . 

1/2 

The PHSD approach was applied to nucleus-nucleus collisions from sJ N ~ 5 to 
200 GeV in Refs. ifTOlfTTl in order to explore the space-time regions of partonic mat- 
ter. It was found that even central collisions at the top-SPS energy of ^/snn =17.3 
GeV show a large fraction of nonpartonic, i.e. hadronic or string-like matter, which 
can be viewed as a hadronic corona lfl6l . This finding implies that neither hadronic 
nor only partonic models can be employed to extract physical conclusions in com- 
paring model results with data. 



3 Application to nucleus-nucleus collisions 

In this Section we employ the PHSD approach to nucleus-nucleus collisions at mod- 
erate relativistic energies. It is of interest, how the PHSD approach compares to the 
HSD [14] model (without explicit partonic degrees-of-freedom) as well as to ex- 
perimental data. In Fig. [3] we show the transverse mass spectra of n~, K + and K 
mesons for 7% central Pb+Pb collisions at 40 and 80 A-GeV and 5% central col- 
lisions at 158 A-GeV in comparison to the data of the NA49 Collaboration lTT7l . 
Here the slope of the n~ spectra is only slightly enhanced in PHSD relative to HSD 
which demonstrates that the pion transverse motion shows no sizeable sensitivity to 
the partonic phase. However, the transverse mass spectra are substantially hard- 
ened with respect to the HSD calculations at all bombarding energies - i.e. PHSD is 
more in line with the data - and thus suggests that partonic effects are better visible 
in the strangeness-degrees of freedom. 
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Fig. 3 The Tl~ , K + and K~ transverse mass spectra for central Pb+Pb collisions at 40, 80 and 
158 A-GeV from PHSD (thick solid lines) in comparison to the distributions from HSD (thin solid 
lines) and the experimental data from the NA49 Collaboration |17| . 
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Fig. 4 The n~ , K + and K~ transverse 
mass spectra for 5% central Au+Au 
collisions at </s = 200 GeV from 
PHSD (thick solid lines) in compar- 
ison to the distributions from HSD 
(thin solid lines) and the experimen- 
tal data from the BRAHMS, PHENIX 
and STAR Collaborations fT8l|T9][20l 
at midrapidity. 



The PHSD calculations for RHIC energies show a very similar trend - the inverse 
slope increases by including the partonic phase - cf. Fig. [4] where we show the 
transverse mass spectra of %~ , K + and K mesons for 5% central Au+Au collisions 
at y/s = 200 GeV in comparison to the data of the RHIC Collaborations Ill81ll91l20ll . 

The hardening of the kaon spectra can be traced back to parton-parton scattering 
as well as a larger collective acceleration of the partons in the transverse direction 
due to the presence of repulsive vector fields for the partons. The enhancement of 
the spectral slope for kaons and antikaons in PHSD due to collective partonic flow 
shows up much clearer for the kaons due to their significantly larger mass (relative to 
pions). We recall that in Refs. (2TJ the underestimation of the K^ 1 slope by HSD (and 
also UrQMD) had been suggested to be a signature for missing partonic degrees of 
freedom; the present PHSD calculations support this early suggestion. 

The strange antibaryon sector is of further interest since here the HSD calcula- 
tions have always underestimated the yield \22\. Our detailed studies in Ref. iflOl 
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show that the HSD and PHSD calculations both give a reasonable description of the 
A+E° yield of the NA49 Collaboration l23l : both models underestimate the NA57 
data |24] by about 30%. An even larger discrepancy in the data from the NA49 
and NA57 Collaborations is seen for (A + E°)/N woum i; here the PHSD calculations 
give results which are in between the NA49 data and the NA57 data whereas HSD 
underestimates the (A + E°) midrapidity yield at all centralities. 

The latter result suggests that the partonic phase does not show up explicitly in an 
enhanced production of strangeness (or in particular strange mesons and baryons) 
but leads to a different redistribution of antistrange quarks between mesons and 
antibaryons. In fact, as demonstrated in Ref. ifTOl . we find no sizeable differences 
in the double strange baryons from HSD and PHSD - in a good agreement with the 
NA49 data - but observe a large enhancement in the double strange antibaryons for 
PHSD relative to HSD. 

The anisotropy in the azimuthal angle y/ is usually characterized by the even 
order Fourier coefficients v„ = (exp( — ffo>))), n = 2,4, since for a smooth 
angular profile the odd harmonics become equal to zero. As noted above, ffo> is 
the azimuth of the reaction plane and the brackets denote averaging over particles 
and events. In particular, for the widely used second order coefficient, denoted as an 
elliptic flow, we have 



where p x and p y are the x and y components of the particle momenta. This coefficient 
can be considered as a function of centrality, pseudo-rapidity 77 and/or transverse 
momentum pj. We note that the reaction plane in PHSD is given by the ix — z) 
plane with the z-axis in the beam direction. 

In Fig.[5]the experimental vi excitation function in the transient energy range is 
compared to the results from the PHSD calculations H301 : HSD model results are 
given as well for reference. We note that the centrality selection and acceptance are 
the same for the data and models. 

We recall that the HSD model has been very successful in describing heavy-ion 
spectra and rapidity distributions from SIS to SPS energies. A detailed compari- 
son of HSD results with respect to a large experimental data set was reported in 
Refs. 127] ED ED for central Au+Au (Pb+Pb) collisions from SIS to top SPS en- 
ergies. Indeed, as shown in Fig. [5] (dashed lines), HSD is in good agreement with 
experiment for both data sets at the lower edge (^/saw ~10 GeV) but predicts an 
approximately energy-independent flow V2 at larger energies and, therefore, does 
not match the experimental observations. This behavior is in quite close agreement 
with another independent hadronic model, the UrQMD (Ultra relativistic Quantum 
Molecular Dynamics) J29) (cf. with Ref. 1231 ). 

From the above comparison one may conclude that the rise of v>2 with bombard- 
ing energy is not due to hadronic interactions and models with partonic d.o.f. have 
to be addressed. Indeed, the PHSD approach incorporates the parton medium ef- 




(3) 
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Fig. 5 (Left: )Average elliptic flow v% of charged particles at midrapidity for two centrality se- 
lections calculated within the PHSD (solid curves) and HSD (dashed curves). The vi STAR data 
compilation for minimal bias collisions are taken from 1251 (stars) and the preliminary PHENIX 
data 1 26 1 are plotted by filled circles. 

Fig. 6 (Right:) The evolution of the parton fraction of the total energy density at mid- 
pseudorapidity for different collision energies. 



fects in line with a 1QCD equation-of-state, as discussed above, and also includes a 
dynamic hadronization scheme based on covariant transition rates. It is seen from 
Fig.[5]that PHSD performs better: The elliptic flow V2 from PHSD (solid curve) is 
fairly in line with the data from the STAR and PHENIX collaborations and clearly 
shows the growth of V2 with the bombarding energy l30l . 

The \>2 increase is clarified in Fig. [6] where the partonic fraction of the energy 
density at mid-pseudorapidity with respect to the total energy density in the same 
pseudorapidity interval is shown. We recall that the repulsive scalar mean-field po- 
tential U s (p s ) for partons in the PHSD model leads to an increase of the flow V2 
as compared to that for HSD or PHSD calculations without partonic mean fields. 
As follows from Fig. [6] the energy fraction of the partons substantially grows with 
increasing bombarding energy while the duration of the partonic phase is roughly 
the same. 

The V2 coefficient measures the response of the heated and compressed matter 
to the spatial deformation in the overlap region of colliding nuclei, which is usually 
quantified by the eccentricity £2 = < y 2 — x 2 > / < x 2 + y 2 > . Since the flow response 
(V2) is proportional to the driving force (£2), the ratio i'2/£2 is used to compare 
different impact parameters and nuclei. 

A remarkable property - universal scaling - has been proposed in Ref. QTl (see 
Fig. [7). It appears that V2/£2 plotted versus {\/S)dN c h/dy falls on a 'universal' 
curve, which links very different regimes, ranging from AGS to RHIC energies. 
Here S = %\/ < x 2 >< y 2 > is the overlap area of the collision system and dN c iJdy 
is the rapidity density of charged particles. 
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• 200 GeV STAR preliminary 

■ 62 GeV STAR preliminary 

A 17 GeV NA49 

T 9 GeV NA49 

+ 5 GeV E877 : 




Fig. 7 Scaling of vi/ej vs 
(l/S)(dN ch /dy). The PHSD 
results |32| are given by lines 
with open symbols. Predictions 
of ideal boost-invariant hydro- 
dynamics are shown in the top 
panel (from 1331 ) and explained 
in the text. Our PHSD results are 
presented in the bottom panel. 
The experimental data points for 
Au+Au collisions at 200 GeV 
(circles) and 62 GeV (squares) 
are from Refs. I33II34I. 



10 15 20 
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As seen fromFig.[7](lowerpannel) the universal scaling of V2/E2 versus (1 /S)dN c i J /dy 
is approximately reproduced by PHSD (see Ref. l32l for the details). This feature is 
not reproduced by hadronic transport models (such as HSD and UrQMD) and meets 
(severe) problems in the various hydrodynamic descriptions as demonstrated in the 
upper pannel of Fig. [7] for a pure hadronic equation of state ('EoS H') as well as 
with a QGP phase transition ('EoS Q'). 

Thus, the experimentally observed scaling in Fig.|7]puts very strong constraints 
on the initial microscopic properties (entropy density, mean free path, etc.), as well 
as the global longitudinal structure l35l . 

Work supported in part by the HIC for FAIR framework of the LOEWE program 
and by DFG. 
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